home *** CD-ROM | disk | FTP | other *** search
/ Celestin Apprentice 2 / Apprentice-Release2.iso / Tools / Languages / RLaB 1.15c / examples / collatz.r < prev    next >
Encoding:
Text File  |  1994-02-21  |  457 b   |  35 lines  |  [TEXT/RLAB]
  1. //
  2. // The 3n + 1 problem
  3. //
  4.  
  5. collatz = function(start)
  6. {
  7.   local(c, n);
  8.   c = n = start;
  9.  
  10.   while(n > 1)
  11.   {
  12.     if(mod(n,2) == 0) {
  13.       n = n/2;
  14.     else
  15.       n = 3*n + 1;
  16.     }
  17.     c = [c, n];
  18.   }
  19.   return c;
  20. };
  21.  
  22. //
  23. // Try it out
  24. //
  25.  c = collatz(100);
  26.  
  27.  plgrid ();
  28.  ptitle ( "RLaB Collatz Example" );
  29.  xlabel ( "Independent Variable" );
  30.  ylabel ( "Dependent Variable" );
  31.  plot( [1:c.nc; c]' );
  32.  pause ();
  33.  plgrid ( ,"bcgnstlv");
  34.  plot( [1:c.nc; c]' );
  35.